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Abstract. The present review is intended to revisit the advances and debates in the 
comprehension of the mechanisms of subcritical crack propagation in sihcate glasses 
almost a century after its initial developments. Glass has inspired the initial insights of 
Griffith into the origin of brittleness and the ensuing development of modern fracture 
mechanics. Yet, through the decades the real nature of the fundamental mechanisms 
of crack propagation in glass has escaped a clear comprehension which could gather 
general agreement on subtle problems such as the role of plasticity, the role of the 
glass composition, the environmental condition at the crack tip and its relation to the 
complex mechanisms of corrosion and leaching. The different processes are analysed 
here with a special focus on their relevant space and time scales in order to question 
their domain of action and their contribution in both the kinetic laws and the energetic 
aspects. 
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1. Introduction 

Silicate glasses are highly homogeneous and isotropic materials at scales larger than a few 
nanometres, thus making continuum descriptions well suited down to submicrometric 
scales. This strong degree of homogeneity is at the origin of the best glass properties, 
such as the excellent transparency, but it is also partly the cause of a major limitation: 
brittleness! Glass products break very easily when subjected to even weak impacts 
with hard materials (Mould 1967), and also the static stress that can be sustained by 
a glass product for a convenient time without fracturing is quite low (Preston 1942, Gy 
2003). The extreme homogeneity of glass along with the elevated strength of its covalent 
cohesive bonds cause a huge stress concentration in a nanometric region around small 
mechanical defects. The local evolution of damage is thus very efficient, involving very 
little bulk dissipation and relatively efficient conversion of the mechanical energy into 
the surface energy required to create new fracture surfaces. 

The earhest work of Griffith (1921) has identified the origin of the low strength of 
glass to the unavoidable presence of small flaws in glass objects. These flaws are generally 
interpreted as submicrometric surface cracks that escape notice under optical inspection. 
Worse than this, these stress concentrators are responsible for a local acceleration of the 
corrosive influence of the environment on the glass. The damage will then generally 
spread in time even under moderate global stresses, leading finally to a delayed failure, 
especially in the presence of a moist atmosphere. This phenomenon, named static 
fatigue, can be interpreted as the combination of an initial stage of differential stress- 
enhanced corrosion (Charles et HiUig 1962), leading to the progressive sharpening of 
the flaws until the crack tip radius reaches a molecular (nanometric) dimension, and 
a second stage of stationary slow crack growth (Wiederhorn and Bolz 1970) with the 
stress-corrosion reactions concentrated at the molecularly sharp crack tip. 

The development of fracture mechanics has allowed a rationalisation of the analysis 
of the fracture resistance of glass by describing the conditions of propagation of a single 
major fracture in a glass sample. The earliest introduction by Griffith of the strain 
energy release rate G, and then of the stress intensity factor K (Irwin and Washington 
1957) have made it possible to determine a more fundamental fracture criterion and 
more robustly defined fracture parameters for crack propagation in the form of critical 
values Gc or Kc- Yet, in most practical cases Gc can be larger than twice the equilibrium 
surface energy of the material, implying the occurrence of other subsidiary dissipation 
mechanisms in the fracture process. 

The brittleness of glass can thus be read as the elevated efficiency of conversion of 
elastic energy into surface energy, suggesting a relatively weak contribution of plastic 
deformation and other bulk dissipation mechanisms, which would involve at most a 
region of few nanometres around the crack tip. The values of Kc for most glasses are of 
the order of 1 MPam^/^, which is rather low, especially when compared to the 50 fold 
larger values in metals where the energy dissipation by bulk ductility is dominant in the 
fracture energy. 
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In the fracture mechanics framework the subcritical propagation by stress-corrosion 
can be well rationalised and represented as a v{K) or G{v) relation for K < 
(Wiederhorn and Bolz 1970). This was shown to generally involve three propagation 
regions (cf. Fig. Hj). Region I corresponds to the cited stress-corrosion regime, where 
the crack propagates due to a stress-enhanced corrosion reaction, which is strongly 
dependent on the environment. In region II the propagation velocity is limited by 
the transport kinetics of environmental corrosive molecules to the crack tip. In region 
III, where the stresses are strong enough to induce the bond breaking in the absence 
of corrosion reactions, the velocity rises very sharply with K, eventually reaching the 
critical propagation at Kc. For some glass compositions, a threshold behaviour appears 
for K = Kq leading to what is generally called region (Michalske 1977). 

Phenomenological equations such as Wiederhorn's (cf. section ISTTI) can describe the 
dependence of crack velocity on stress and on environmental parameters for most typical 
glasses and are compatible with a thorough consistent modelling based on the sharp- 
crack atomic-bonding paradigm (Lawn 1993). Yet the detailed nature of the stress- 
corrosion mechanisms that occur at the crack tip have been debated for decades (Marsh 
1964b; Maugis 1985; Gehrke et al 1991; Tomozawa 1996), and a general disagreement 
can be found on the relevance of several accessory phenomena that may participate in 
the stress-corrosion mechanisms at different stages of the process. Stress-corrosion can 
involve a complex interplay between the diffusion of reactive molecules (mainly water) 
into the crack cavity and into the glass network, the corrosion (or dissolution) of the 
network itself, and the migration of weakly bonded alkali ions under chemical or stress 
gradient (Gehrke et al 1991; Bunker 1994). All these phenomena are typically very 
slow under ambient conditions in the unstressed material, but they can significantly 
accelerate in the highly stressed neighbourhood of the crack tip, depending especially 
on the nature of the environment and of its confinement. Their time scales will be 
progressively accelerated in a series of smaller and smaller shells surrounding the crack 
tip (cf. Fig. [1]), defined by the stress being larger than some critical value. However, 
since the crack tip is propagating, the different phenomena will be able or not to spread 
in each shell depending on a competition between their time scale in each shell and 
the crack propagation velocity. Since the size of the activated shells can vary from 
microns to nanometres depending on the interplay between space and time scales, the 
modelling of these phenomena (which will be reviewed in more detail in next sections) 
requires nanometre scale resolved investigation techniques and questioning about the 
relevant physical laws at the nanoscale. Such studies should be very promising to solve 
the debates by direct observation and to relate the phenomenological parameters to the 
specific composition and structure of glasses. 

The development of several nanoscale investigation techniques in the 80's has 
led to significant insights into the comprehension of the different mechanisms (cf. 
sections 14.31 to W^ . In particular SEM, TEM and AFM measurements of both post- 
mortem crack surfaces and direct in-situ observations of the crack tip neighbourhood 
have permitted the investigation of the space and time scales of occurrence of these 
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phenomena under specific conditions. The progressive increase in resolution down to 
the micron scale of several structural and compositional analysis techniques (such as 
Raman, IR, Brillouin spectroscopies, X-Ray, electron and neutron scattering, NMR, 
XPS, SIMS) have also permitted investigations of the alterations of the bulk glass near 
the crack tip or at freshly fractured surfaces. These techniques, combined with the 
increasing power of molecular and quantum dynamics simulations permit great insights 
into the understanding of the combination between mechanical and chemical processes 
acting at the crack tip. 

Section [2] will first focus on some relevant concepts of fracture mechanics. In section 
[3] the most solidly established features of the slow crack propagation kinetics in glasses 
will be described and interpreted in the framework of the classic sharp-crack atomic- 
bonding paradigm proposed in the '70s. The section|4]of the review will present a deeper 
discussion of the relations between the crack propagation and the chemical mechanisms 
of stress-corrosion, along with a critical analysis of the points which are still controversial 
and of the recent experimental evidences and efforts which have been done to make 
them clearer. A concluding section will discuss the perspectives and promises of the 
development of this interesting research field. 

2. Space, time and energy scales of fracture mechanics 

The world of fracture mechanics is extremely wide and rich of theories and variants that 
account for the many peculiar manifestations in different materials. It is difficult to 
figure out the 'one theory' that can encompass different phenomena such as atomic bond 
breaking, dislocation motion, viscoelastic-plastic flow, damage spreading, corrosion, ion 
exchange, cavitation, fingering and a whole plethora of other fancy effects. Different 
application communities generally end up with following different paths, along with 
different terms, variables and philosophies. Luckily enough, the classic Linear Elastic 
Fracture Mechanics (LEFM) was mainly developed to tackle strength problems of brittle 
materials, and glass has often been a privileged material for its understanding and 
application (cf. Lawn, 1993). In the present section, I will revisit the main relevant 
concepts and theories of fracture mechanics and propose a synthetic point of view that 
will guide the analysis in the present review by focusing on the space, time and energy 
scales of the different physico-chemical phenomena occurring during fracture. 

Although the ideal basic mechanism for fracture is the conversion of mechanical 
energy into surface energy (Griffith 1921), fracture propagation is generally complicated 
by the activation of subsidiary dissipation mechanisms that generally act in a bulk 
neighbourhood of the crack tip called the process zone (Irwin 1960), but that can also 
extend to the whole material volume or act between the open crack walls as illustrated in 
the schematic model in Fig. [H The energy dissipated by each of these mechanisms has to 
be taken into account to estimate the energy per unit area that is necessary to propagate 
the fracture. To progressively account for nonlinear effects, irreversibility and dissipative 
processes, fracture mechanics has undergone several profound developments, which 
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however were managed with appreciable devotion in order to preserve the formahsm 
of LEFM (Orowan 1955; Irwin 1960; Barenblatt 1962; Rice 1968; Rice 1978; Lawn 
1993). This is based on one hand on the use of surface energy balance related variables 
such as the strain energy release rate G and on the other hand on variables describing 
the fundamental inverse square root dependence of the stress field on the crack tip 
distance such as the stress intensity factor K (Irwin and Washington 1957)0. Moreover, 
the two kinds of variable remain in general related by an equivalence in the typical form 
G = K^/E' {E' being an equivalent Young's modulus) and the criteria for equilibrium 
and crack propagation can still generally be expressed by equations of the kind G > Gc 
or K > Kc by only changing the exact definition of the variables used (Lawn 1993). 

As concerns the description of the kinetics of fracture, we can distinguish two main 
families of approaches: the ones where the dissipated energy dominates the surface 
energy and the ones where the crack tip debonding processes play the dominant role. 
These approaches involve the use of different variables and diagrams, but also reflect 
different philosophies. In order to provide a common frame it is useful to analyse the 
mechanical scheme in Fig. [TJ In a typical fracture propagation experiment, the external 
load can be represented by a given force F, displacement X or traction velocity V 
applied to a loading point. In the most general case, the test sample can be subdivided 
into a series of regions that progressively transfer the energy towards the crack tip, 
such as the body, the process zone (one or a series of shells depending on the number 
of potential dissipative mechanisms), a local near-tip elastic enclave, and the region 
inside the crack cavity. The mechanical role of each of these regions is schematically 
represented in the block diagram in Fig. [21 where the nature and spatial extension of each 
block can change between different materials and depend on the time scales of loading 
and crack propagation, each block possibly being elastic, viscoelastic, viscoplastic, etc 
(the symbols used in the graph are just evocative!). In order for K to be definable and 
the LEFM to apply, the body block should be purely elastic and the process zone size 
should be small with respect to the sample and fracture length sizes. According to the 
theory for brittle materials (Lawn 1983), it is possible to define a small enclave zone near 
the crack tip that behaves in an essentially elastic way and that allows the definition of 
the local values of G* and K* that are relevant for describing crack tip processes. 

In the first approach (mainly formalised by Maugis (1985) to describe fracture in 
polymers and peeling of adhesives) fracture propagation in a given material is generally 
characterised by G{v) curves, representing the strain energy release rate as a function 
of the crack velocity (Fig. [3]). This makes it particularly easy to consider the energy 
balance of fracture propagation in relation to the time scales of deformation solicited 
by different crack propagation velocities. In this formalism, the low velocity limit can 
be identified with the 'equilibrium' surface energy Gq = 270. Then the fracture energy 
increases with velocity due to enhanced activation of concomitant irreversible dissipative 
processes, which determine a stationary propagation curve in the form G{v). For the case 

I The present review is limited to mode I tensile fracture for simplicity, the variable K will stand for 
Ki 
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Figure 1. Graphic representa- 
tion of the mechanical elements 
for fracture mechanics. 



Figure 2. Symbolic block 
representation of the different 
mechanical elements. 



of viscoelastic materials, the slow fracture region is generally in the form of a power law 
of the kind G{y) = Go{l + Atv""), where At includes the temperature dependence. The 
increasing portion of the G{v) curve eventually ends at some critical velocity Vc, where 
the dissipative process starts to be less and less effective, since the kinetic scale becomes 
faster than the viscoelastic relaxation. For higher velocity the fracture energy increases 
again due to the manifestation of dynamic effects (Freund 1990) eventually leading to 
a divergence when the crack velocity approaches the sound propagation velocitj|§|. The 
presence of a maximum in the G{v) curve leads to an intrinsic mechanical instability in 
the fracture propagation dynamics for a critical value Gc that should not be confused 
with the threshold surface energy Go = 270, although it is generally proportional to 
it. The range of velocities corresponding to a negative slope in the G{v) curve does 
not allow stationary crack propagation and generally leads to stick-slip dynamics (cf. 
Ciccotti et al 2004). The instability at Gc, i.e. the jump from slow to rapid fracture for a 
constant applied force, is of particular interest since it is marked by a sharp reduction of 
the amount of energy dissipated by fracture propagation, letting all the excess furnished 
energy to be converted to kinetic energy with the typical emission of acoustic bursts 
(Barquins and Ciccotti 1997). 

In the second approach (mainly formalised in the work of Wiederhorn et al (1982) 
on glass and Lawn (1993) for brittle materials) the focus is on the kinetics of a crack 
tip debonding process as a function of the applied stresses. The fracture propagation 
curves are generally represented by v{K) curves, where the crack propagation velocity 

§ The observed limit propagation velocity is generally only a fraction of the sound velocity due to the 
occurrence of other kinds of dynamical instabilities (cf. the review of Fineberg and Marder 1999). 
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Figure 3. Typical G{v) diagram 
for crack propagation in poly- 
mers and adhesives (w = Gq). 



Figure 4. Typical v(K) or v{G) 
diagram for subcritical crack 
propagation in glass. 



is considered as a variable dependent on the stress intensity factor (Fig. H]). The 
global energy balance is somehow in a subsidiary position, and the fracture kinetics is 
determined by the rate of thermal activation of a generic discrete bond breaking reaction 
(Thomson et al 1971) whose activation energy is reduced by the application of stresses, 
leading to an Arrhenius like exponential equation v = foexp[— (Ai^^ — bK)/kBT], that 
may also be expressed in terms of G (cf. section [3?T]) . In this interpretation, the critical 
value Kc for the transition to dynamic propagation should correspond to the end of the 
lattice trapped regime, while the equilibrium surface energy 270 should correspond to 
a propagation threshold. However, this basic modelling can only be argued for a single 
debonding mechanism. The inert fracture (cf. section 13. 3p should thus be modelled 
separately from chemical assisted propagation (cf. sections 13.11 and 13.41) . We should 
however retain the idea that the position of the critical value Kc is related to the 
discreteness and heterogeneity of the material rather than to a surface energy, since it 
represents the moment where the energy barriers for individual bond breaking tend to 
zero. The surface energies should rather be compared with the propagation thresholds 
for each mechanism. In materials which are heterogeneous on a large scale, this approach 
can be reformulated in terms of a larger basic grain size (Broberg 2003). 

The two approaches can be usefully combined in brittle materials presenting a non 
negligible process zone, but where a small near-tip elastic enclave can be defined. To 
do that we should first assume that the local values of G* and K* in the enclave are 
effective in controlling the crack tip propagation kinetics according to a relation v{K*) 
(Lawn 1983). In the absence of forces inside the crack cavity, the energy dissipated in 
the process zone can be evaluated by the difference G — G* (applying the J integral 
formalism of Rice 1968) and will cause both G* and K* to be smaller than their 
macroscopic counterpart G and K. The effect of energy dissipation will thus translate 
into a shielding (or local relaxation) of the crack tip stresses. The presence of generalised 
forces in a limited portion of the near-tip crack cavity will also act as a shielding factor 
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at the crack tip, and it may or may not involve energy dissipation depending on the 
nature of the interaction. 

When describing the macroscopic fracture properties of a material, it is the relation 
of the external variables G, K which is studied as a function of v. Even if assuming a 
uniquely defined kinetic law v{G* ,K*), these local variables remain hidden, and the 
presence of a process zone will cause the non-uniqueness of the macroscopic fracture 
laws f (G, K). This induces a time- or crack-length dependent behaviour (R-curve) that 
may lead into controversial interpretation (Lawn 1993). These types of modelling will 
be shown later to be particularly effective in rationalising the propagation threshold 
behaviour (cf. section [331). 

As a concluding remark to this section, it is generally useful to identify the different 
space, time and energy scales of each symbolic block (cf. Fig. [2]) in order to estimate 
if the different dissipation mechanisms will have time to activate and how they will 
influence the energetic balance or the fracture propagation laws. Several spatial scales 
are present in the problem, such as the sample lateral dimensions and thickness, the 
loading displacement, the crack length, the basic grain size of the material (molecular 
rings for glass), the size of the process zone (or of each shell of a multiple process zone), 
the enclave size, the length of the one or several cohesive zones in the crack cavity. 
Typical time scales are given by the loading or displacement rate, the crack velocity, the 
stimulation time of each equivalent block and its characteristic relaxation time, the rate 
of thermal activation of different mechanisms, the rate of transport of different relevant 
chemical species in the bulk or in the crack cavity. We will try to keep these space and 
times scales in mind while proceeding in the description of the crack propagation kinetics 
in glasses and especially when trying to get a deeper insight into the stress-corrosion 
mechanisms. 

3. Crack propagation kinetics 

A huge variety of possible oxide glasses may be created by melting together variable 
amounts of oxides, within specific compositional ranges that are strongly influenced 
by the quench speed (cf. Zarzycki 1991). The details of the chemical reactions and 
the competition between different propagation mechanisms will be altered, but many 
typical features of crack propagation remain substantially similar. The most influential 
difference is between glasses essentially made by network-former oxides {Si02, Ge02, 
B2O3, P2O5, etc.) and glasses containing significant amounts of modifier oxides {Na20, 
K2O, CaO, etc.). For that reason in the following discussion silica glass will mainly be 
described as a representative of the first class and a typical soda-lime glass will be used 
as paradigm to represent the second class of glasses. 

The subcritical fracture propagation properties of glasses are generally studied using 
samples that allow the stable slow propagation of a single fracture such as the Double 
Torsion (DT, Evans 1972), the Dual Cantilever Beam (DCB, Srawley and Gross 1967) 
or the Double Cleavage Drilled Compression (DCDC, Janssen 1974). The measurement 
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of the crack propagation velocity as a function of the apphed stress intensity factor 
generally presents three (or four) characteristic regions, Fig. |H 



3.1. Region I: stress- corrosion regime 

In region I, corresponding to the stress-corrosion regime, the crack propagation velocity 
is a strongly increasing function of the stress intensity factor K; it has an almost linear 
dependence on humidity in moist air (Fig. [5]) and it increases with temperature (Fig. 
|6]). All these dependencies can generally be fitted by the Wiederhorn (1967) equation: 

V =VQe^^[aK) = A\——\ exp I — I (1) 

where PH2O is the partial pressure of the vapour phase in the atmosphere, po is the 
total atmospheric pressure, and R the gas constant; A, m, AEa and b are four 
adjustable parameters that take into account the dependence on the glass composition 
(cf. interpretation below). 




Figure 5. Effect of humidity on 
the crack propagation in soda- 
hme glass (from Wiederhorn 
1967). 



Figure 6. Effect of temperature 
on the crack propagation in soda- 
lime glass (from Wiederhorn and 
Bolz 1970). 



However, due to the extremely strong dependence of velocity on the stress intensity 
factor K and to the limited range of variation of K, the data in region I may also be 
fitted by a power law expression like in the Maugis (1985) approach: 



= VoiK/Ko) 



(2) 
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the exponent n being generally between 12 and 50 for silicate glasses (Evans and 
Wiederhorn 1974). This form is by the way particularly useful (and indeed used) for 
developing simple analytical predictions of the life times in static and dynamic fatigue 
tests (Davidge et al 1973). The difficulty in determining the most adequate relation 
results in a strong uncertainty in the predictions over very long periods. 

For the same reasons, the region I data can also be fitted by a different exponential 
equation proposed by PoUet and Burns (1977) and applied to glasses by Lawn (1993): 

V = t;oexp(«C;) = Aexp ^-^j exp I -j-^ I (3) 

where the energetic balance has a more central role, as will be discussed later concerning 
the threshold at Go- 




Figure 7. Basic mectianisni 
of tiie stress-corrosion reaction 
(from Michalske and Freiman 
1983). 



Figure 8. Effect of glass compo- 
sition on the crack propagation 
in soda-lime glass (from Wieder- 
horn and Bolz 1970). 



The Wiederhorn equation ([T]) has been the most used in the description of region I in 
glasses since it provides a more direct interpretation of the chemical reactions between a 
glass of a given composition and the environment at the crack tip. In this framework, the 
basic mechanism of the stress-corrosion reaction was associated with the stress-enhanced 
thermal activation of a dissociative hydrolysis reaction represented in Fig. [7] for the case 
of silica glass (Michalske and Freiman 1983). Michalske and Bunker (1984) supported 
the positive effect of stress on the reaction rate by a molecular-orbital simulation of the 
reaction on distorted siloxane bonds. They later provided experimental support for their 
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theory by measuring the rate of hydrolysis of a series of strained sihcate ring structures 
(Michalske and Bunker 1993). For other glass compositions (Fig. [8]), several other 
reactions are possible, but they share the similar feature of being stress-enhanced and 
thermally activated, thus explaining the exponential (Arrhenius like) term in equation 
([I]), where A£^a can be interpreted as an activation enthalpy in the absence of stress 
and h can be expressed in terms of an activation volume (of molecular dimensions) by 
relating the stress intensity factor to the crack tip stress (Wiederhorn and Bolz 1970). 
According to the chemical reaction rate theory (Glasstone et al 1941), the preexponential 
term represents the activity of the reactants in the moist atmosphere and the exponent 
m (generally close to 1) is related to the nearly first order of the stress-corrosion reaction 
(Wiederhorn et al 1974). 

For propagation in liquid environment, the term p/po is substituted by the activity 
of water molecules which is strongly dependent on the composition of the liquid 
environment. Notably, it is strongly affected by the pH and ionic concentration of 
some species related to the glass composition (Wiederhorn and Johnson 1973). Apart 
from water, other molecules that are effective for stress-corrosion such as ammonia, 
hydrazine and formamide have in common (a) the capability of donating both a proton 
and an electron at the two ends of a molecule and (b) a molecular diameter inferior 
to 0.5 nm (Michalske and Bunker 1987). The first property is necessary to present a 
concerted reaction of adsorption and scission of the siloxane bridges of the glass network 
(Fig. [71 bottom). The small diameter is necessary for the reactive molecules to be able 
to reach the strained bonds at the crack tip without steric problems. The effect of other 
non reactive liquids is mainly that of reducing the effective concentration of water to 
traces and thus of substantially suppressing the stress-corrosion propagation (Freiman 
1975) to lead to a propagation similar to that in dry nitrogen or in vacuum. 

3.2. Region II: transport limited propagation 

In region II the crack propagation velocity is substantially independent of the stress 
intensity factor (Fig. [5]), since the reaction rate is limited by the transport kinetics of 
reactive molecules towards the crack tip. The plateau value for propagation in moist air 
is increased when increasing humidity in qualitative accord with the proposed equation 
(Wiederhorn 1967): 

= V0PH20DH2O (4) 

where DH2O is the diffusion coefficient of water molecules in air. However, a finer 
analysis reveals that several mechanisms can act in limiting the transport of the reactive 
molecules as a function of the environment (total pressure, humidity, temperature, liquid 
or vapour form of the water molecules) and of the progressive confinement experienced 
when approaching the crack tip (Lawn 1993). The diffusion of water molecules in a 
vapour or solute will be first reduced by the decrease of the average free path, and then 
possibly by the activated adsorption on the crack walls when the confinement becomes 



CONTENTS 



13 



of the same order as the molecular size. While the first mechanism is generally accepted 
for glasses (Michalske and Bunker 1987), the second one is likely to be determinant in 
the cleavage of mica under humidity (Wan, Aimard et al 1990). 

Region II is very limited for samples plunged in liquid water, and the extension of 
the region II in glass was shown to be reduced at elevated humidity (Wan, Lathabai 
et al 1990) suggesting the presence of capillary condensation in the crack cavity. In 
liquid solutions, the position of region II depends on the water concentration, and a 
region II plateau remains generally visible in a very low velocity region for inert liquid 
or gas environments due to trace amounts of water as discussed above (Freiman 1975). 
Water molecules may also penetrate the strained glass network in a molecular form 
(Doremus 2002) and react with the network in a more internal location near the crack 
tip (Marsh 1964b; Tomozawa 1996). In this case the bulk diffusion rate could also act 
as a limiting factor in region II. However, the kinetics of most of these phenomena is 
ruled by the thermal activation of reactions with a stress dependent energy barrier, and 
is thus expected to depend on temperature and stress in a similar manner (cf. section 
14. ip . making the discrimination difficult if v{K) curves are only measured. 

3. 3. Region III: inert propagation and failure 

When reaching region III, the crack propagation velocity starts increasing again with K 
(Fig. [5]). The dependence is so steep that the measurements require special techniques, 
such as high speed cameras or a modulation of the fracture surface generated by sound 
waves (Gehrke et al 1991). The position of region III is substantially independent 
of the local environment (except for a dependence on the dielectric constant as 
shown by Wiederhorn et al 1982), in agreement with the expected suppression of 
the stress-corrosion reactions. The propagation in vacuum or in inert liquid or gas 
environment extends the region III domain down to very low velocities (Wiederhorn et al 
1974). These delicate measurements showed that the v{K) curves are still temperature 
dependent and that they can be fitted with an Arrhenius relation with a more elevated 
activation volume. 

Although the velocity is rapidly climbing to large values, the propagation is still of a 
subcritical type, i.e. the time scale for crack propagation is still determined by the time 
required for thermal fluctuations to overcome the weaker and weaker energy barriers. 
According to lattice trapping models (Thomson et al 1971), the critical stress intensity 
factor Kc represents the end of the region III, i.e. the moment where no more energy 
barriers exist and the mechanical energy is converted into kinetic energy in an unstable 
manner, the time scale for propagation being now determined by dynamic equations. 
However, due to the elevated slope of region HI in glasses and to the high values of the 
crack velocity, it is hardly possible to establish experimentally the end point of region 
III and the value of Kc is generally identified with the position of region HI itself. From 
a practical point of view, when loading a sample in tension, the sample will generally 
fail abruptly when approaching Kf., and the measured value of K at failure is called the 
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fracture toughness. On the other hand, in stable fracture configurations [dK/dc < 0, 
c being the crack length) such under DCDC conditions, and in the absence of stress- 
corrosion (for example in vacuum or in nitrogen) the crack will simply arrest for K = Kc 
and further propagate only upon an increase of the load. 

The toughness of silicate glasses is generally around 0.8 MPam^(^ within a range 
of ±20% for most typical glass compositions (Vernaz 1978, 1980)[|jl. The relation 
between Kc and glass composition is quite complex since the glass network structure 
is strongly affected by the percent concentration of each oxide. Since the toughness is 
generally related to the Young's modulus, according to the estimation Kc = y/'jE, the 
compositions with higher Young's modulus are generally tougher. Another apparently 
consistent correlation appears between the decrease of toughness and the number of 
non-bridging oxygen atoms, which is determined by the concentration of modifier oxides 
(Rizkalla et al 1992). Moreover, the toughness tends to increase when the size of the 
alkali ion in ternary glasses decreases (Gehrke et al 1991), suggesting the importance 
of ion migration even under the rapid fracture (cf. section 14. 5p . However, a clear 
comprehension is still lacking, and several research efforts are currently underway to 
establish clearer correlations (Kingston and Hand 2000; Kurkjian et al 2003). The 
major reason for this delay is that the measurements of K^. are quite delicate and that 
they are generally obtained with very different techniques, the overall scattering being 
comparable to the strength of the observed concentration effects. 



3.4- Region 0: threshold behaviour, healing and aging 

In some glasses, such as alkali rich glasses, a threshold behaviour occurs (Fig. [9]), i.e. 
fracture velocity rapidly falls to undetectably low values when K decreases toward a 
threshold value K^ (Gehrke et al 1991). This property is of extreme importance in the 
life duration of engineering glass products, yet its origin is very subtle and has long been 
subject to debate. 

According to Griffith theory, a threshold for propagation should exist at G = Go 
where Gq = 270 represents the thermodynamic surface energy which is necessary to 
create the new surfaces through the breaking of the cohesive bonds. In the framework 
of the rate reaction theory (Glasstone et al 1941) this equilibrium condition corresponds 
to the equality of the rates of the two opposing reactions of breaking and recombination. 
The crack propagation velocity in proximity to such an equilibrium threshold can thus 
be approximated by including in Eq. ([3]) a term for the rate of recombination (Lawn, 
1993): 

/ AEa~a{G-Go)\ , / AE, + a{G-Go)\ 
V = .oexp ^ j + .oexp ^ j (5) 

II We should remark that the corresponding fracture energy of glass in inert atmosphere is of the order 
of Gc — 27/ = K"^/ E' ~ 8 J/m^, which is ten times higher than the typical values of the surface tension 
of glasses 70 — 0.5 J/m^ (cf. sect. 14. 3p . 
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Equation ([5]), besides presenting a threshold for propagation at G = Go, also implies 
the backward propagation or healing of the crack for G < Go, in agreement with the 
thermodynamic considerations of Rice (1978)0. 

However, this scenario is not exactly what is happening with glass. For soda-lime 
glass, although the crack propagation stops at a threshold Gg — 0.8 J/m^, no healing 
is observed down to a lower value Gc — 0.15 J/m^ (cf. Fig. [10]), suggesting that the 
crack closing occurs due to a different mechanism from crack opening, for instance 
through the hydrogen bond interaction between hydrolysed crack surfaces (Stavrinidis 
and HoUoway 1983). For silica glass, although no propagation threshold Gg is observed 
down to V = 10~^^ m/s (Muraoka and Abe 1996), crack closure also happens for a value 
of Gc similar to that of soda-lime glass. Moreover, the repropagation kinetics is quite 
different in the two kinds of glasses. While repropagation of the healed crack in silica 
happens as soon as G > Gc, the repropagation threshold G'^ in soda-lime progressively 
increases between Gc and Gg as a function of the time spent in the closed configuration 
(Stavrinidis and Holloway 1983), suggesting partial reformation of the siloxane bridges 
(Michalske and Fuller 1985). 

This phenomenon should not be confused with another kind of aging that can be 
observed in soda-lime glass. When the crack is arrested and held for a given time at a 
value Gh in the interval Gc < Gh < Gg (NB: without healing), the repropagation will 
not start again at Gg, but it will either need the application of a higher value of G, or it 
will occur after a delay time which is maximum when Gh — Gg (Michalske 1977, Gehrke 
et al 1991). 

The threshold behaviour and the delayed fracture after crack arrest in alkali 
containing glasses were shown to be fundamentally related, and can be interpreted as 
a consequence of the ion exchange phenomenon (Gehrke et al 1991; Fett et al 2005a). 
Ion exchange at the crack tip and at crack surfaces can have several effects: (a) the 
build up of compressive stresses at the crack tip, shielding the external stresses, (b) a 
change of the surface energy and/or a dissipation of energy that can alter the energetic 
balance, (c) a change of the local glass composition leading to an increase of the fracture 
resistance. The ion exchange interpretation was supported by the strong dependence 
of the threshold and arrest time for soda-lime in liquids on the pH of the solution and 
the glass composition in a series of binary and ternary glasses (Fig. [9]). Thus, an acidic 
solution (HCl) enhances the ion exchange and causes the threshold to shift to higher 
values of Gg, while a basic solution (NaOH) can make the threshold drop below the 
detection limit. A specific exception is observed in glasses with high alkali content, where 
the corrosion behaviour is so rapid that fracture propagates spontaneously, presenting a 

% We should remark that we can only talk about reversible crack propagation for G « Go, i.e. for 
vanishing propagation velocity, while both forward and backward propagation with finite velocity imply 
a rate of energy dissipation (G — Go)v > 0. 



CONTENTS 



16 




Figure 9. Threshold be- 
haviour in a sodium-alumino- 
siUcate glass in water and in so- 
lutions with different pH (from 
Gehrke et al 1991). 



Figure 10. Sketch of the 
hysteretic healing and reprop- 
agation behaviour in soda-lime 
glass. 



horizontal plateau on the G{v) curve for a low limit velocity around 10~^ m/s, especially 
in presence of an acidic solution (Gehrke et al 1991). 

The phenomenon of ion exchange will be discussed in more detail in sections 14.11 
and 14. 5[ firstly because the effect of local crack tip stresses plays an important role, and 
secondly because the discussion of the difference of behaviour between the propagation 
in water and in moist atmosphere requires a more detailed analysis of the chemical 
processes in confined liquid films. However, some considerations are opportune here. 
The hysteretic behaviour of the crack closure and repropagation in soda-lime glass, 
as well as the non observation of any threshold in silica, cast strong doubts on the 
applicability of the Griffith scenario for reversible healing in glasses. The hydrolysis 
reaction of the siloxane bonds is essentially irreversible at ambient temperature, and 
the crack closure at Gc principally involves the adhesion of the hydrated crack surfaces 
(the latter being fairly reversible for silica). For alkaline glasses, the modifications of the 
crack surfaces are even more drastic due to the irreversible production of ion exchange 
layers. It is therefore not easy to associate the position of the threshold Gg with a glass 
surface energy in the Griffith sense, but its value appears to arise from a competition 
between the kinetics of the formation of the ion exchange layer and the crack propagation 
kinetics. It is however difficult to disregard the similarity between the typical values 
of Gg ~ 0.8 J/m^ and the estimates of the glass surface tension, both being about 
one tenth of the fracture energy in inert environment. This seems to suggest that the 
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complex interaction of glass with water during slow fracture helps the glass surface to 
relax from the elevated non equilibrium energies towards an equilibrated condition. 

4. A deeper insight into the stress-corrosion and damage mechanisms 
4.1. Competition between dissolution, corrosion and leaching 

In order to understand more fully the different components of the chemical contribution 
of water molecules to the corrosion and damage of glass, the three main categories of 
reactions that can occur between glass and an aqueous solution should be reexamined 
(cf. Bunker 1994 for a more detailed review): (1) hydration consists of the penetration 
of molecular water into the glass network as an intact solvent; (2) hydrolysis and 
condensation are the two opposite chemical reactions which break or reform the 
Si — O — Si (or equivalent) network bonds by exchange with a couple of hydroxyl 
groups according to the scheme: 

Si-0-Si + H2O ^ Si- OH + HO - Si 

(3) ion-exchange reactions consist of the replacement of a modifier cation (such as 
sodium) by a hydrogen or hydronium ion: 

Si-0- Na+ + HsO+ ^Si-OH + iVa+ + H2O 

A good account on the solubility of silicate glasses can be found in Her (1979). 
The dissolution reactions are so slow at ambient condition that glasses are considered 
as exceptionally inert materials, and generally need tens of years exposure to moisture 
before showing appreciable surface alteration. However, when increasing temperature 
or stress these reactions can be significantly accelerated, especially under extremely acid 
or basic environmental conditions. Moreover, if the three reactions occur at the same 
time, each reaction will infiuence the kinetics and mechanisms of the other reactions 
(Bunker 1994). 

In pure silica glass only the first two reactions are possible and they cooperate in the 
water difi'usion kinetics. According to the reactive diffusion model developed by Dorcmus 
(well reported in his book of 2002), water mainly diffuses in molecular form with a well 
defined thermally enhanced diffusion coefficient, but it continuously reacts with the 
glass network being thus temporarily immobilised in the SiOH form. As a result, an 
effective water diffusion coefficient can be defined that also depends on concentration. At 
elevated temperature, the water content converges to the equilibrium solubility (mainly 
dependent on the free volume in the glass network) and the hydrolysis reaction is 
equilibrated (mainly towards the dissociated state). However, at lower temperatures 
close to ambient condition, the reactions are not equilibrated and both hydroxyl and 
molecular water arc present in the glass with a very complex speciation, which is often 
subject of debate (Stuke et al 2006; Doremus 2002). 
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In alkali containing glasses such as soda-lime, water diffusion is greatly complicated 
by the combination with the ion exchange reactions (Doremus 1975). Since alkali ions 
partly fill the network voids, ion exchange can be the limiting factor for water diffusion, 
but hydrolysis is also necessary to open up the pathways in the glass network (Bunker 
1994). Lanford et al (1979) have shown that alteration layers are formed on glass 
surfaces in liquid water by selective leaching of alkah ions. The layer width is initially 
proportional to the square root of time (diffusion controlled process) and eventually 
reaches a stationary value due to the equilibration with the slow etching of the glass 
surface. The hydrated layers that are formed on soda-lime glass were shown to be 
sodium-depleted in a proportion that suggests the substitution of sodium ions by the 
hydronium ions H^O^ (rather than the hydrogen ion). This exchange is responsible for 
the development of a stressed surface layer depending on the relative volume difference 
of the exchanged ions. In the case of soda-lime glass the developed stress is compressive 
(Michalske and Bunker 1989, Fett et al 2005b), in agreement with the hypothesis of a 
sodium/hydronium exchange, the volume of hydronium and hydrogen being respectively 
larger and smaller than that of sodium. In an analogue manner the exchange of sodium 
with the larger potassium ions is induced in the chemical tempering treatments to 
enhance the strength of glass products (see Gy 2008 for an extensive review). The 
composition and structure of the leached layer are thus heavily altered and are more 
similar to a hydrated silica glass than to the initial glass. Moreover, under extensive 
leaching, the recondensation reactions can also become important actors in a more 
profound restructuring of the glass network leading to the formation of a more compact 
glass phase, similar to anhydrous silica glass, pervaded by a more permeable porous 
network, often called a gel layer (Bunker 1994). 

While the behaviour and thermodynamics of these corrosion processes have been 
extensively studied as functions of temperature, pH and composition of the liquid in 
contact with the glass surface (cf. the review of Conradt 2008), the behaviour under 
moist atmosphere is often more subtle (Bunker 1994) due to the formation of thin 
nanometric water films on glass surfaces whose detailed chemistry escapes state of the 
art thermodynamics and often entails the formation of complex patterns (Watanabe et 
al 1994). Moreover, although these reactions are known to be strongly infiuenced by the 
local stress state, the accurate study of their behaviour at elevated tensile stresses is still 
lacking and is subject of strong debate. In the next sections we will try to clarify the 
present knowledge about the damage processes near the crack tips, where the extremely 
elevated stress gradients affect smaller and smaller zones, thus making the space and 
time scales of the different mechanisms be highly variable and interdependent. 

4-2. Crack tip blunting 

The initial sharp-crack atomic-bonding paradigm really means that at room temperature 
the hydrolysis of individual crack-tip bonds, which is singularly enhanced by the elevated 
crack tip stress, is the only relevant process for determining the crack tip propagation 
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laws (cf. Fig. [7]). Alternatively, this condition at least should hold once we have identified 
the stress intensity factor K* in the local enclave, thus admitting the possibility that the 
other processes and reactions are effective in a surrounding process zone (Fig. [1]). This 
hypothesis constitutes a singular (degenerate) version of the original Charles and Hilligs 
(1962) theory, in which the stress-enhanced corrosion rate can not act any further in 
sharpening the crack tip since its curvature radius has reached a lower limit, estimated 
to 0.5 nm and corresponding to the radius of the basic siloxane rings of the glass network. 

The assumed condition of propagation at constant tip radius is the basis for the 
interpretation of the Wiederhorn equation ([1]) in terms of the LEFM theory and for the 
uniqueness of the v{K*) curve (Lawn, 1983). However, it can not apply to the description 
of the initial development of blunt flaws into sharp cracks, and we can rise the question 
of whether this evolution would be the only possible destiny for progressively sharpening 
cracks under all conditions. Can the flnal crack tip radius of curvature be a function 
of K at least for certain glass compositions and/or environments? Can the competition 
between the rate of the different processes and the crack tip velocity make the crack 
tip so blunt (or 'effectively' blunt due to material damage) that it may stop under an 
apparently elevated value of K (i.e. larger than that expected from the glass surface 
energy)? The debate on this subject has been especially active concerning the origin 
of the propagation threshold and of the crack closure/aging/repropagation behaviour 
(plasticity issues will be discussed in section H73ll . 

Since the solubility of glass surfaces is also dependent on the local radius of 
curvature (Her 1979), when the external stresses are very low the crack walls should 
be corroded more rapidly than the tip, leading to a progressive tip blunting effect (Ito 
and Tomozawa 1982). Moreover, the presence of a gradient of solubility in a very 
confined environment can lead to a concomitant phenomenon of dissolution from the 
lateral walls and reprecipitation at the crack tip which would even enhance the extent 
of the tip blunting (cf. Fig. [TTl) . as observed in a controversial TEM study of the 
crack tip in silica (Bando et al 1984). If such blunting happens under almost static 
conditions, then the fracture should be substantially reinitiated after reloading, leading 
to an effect of temporary increase of the glass strength, which was proposed as an 
alternative explanation to the occurrence of a time delay before repropagation of aged 
fractures (Han and Tomozawa 1989). 

The difficulty in this debate arises from different causes. On one hand, since all 
these reactions are influenced by stress, temperature and pH, the extent and rapidity 
of occurrence of these phenomena is strongly dependent on the speciflc conditions and 
on both the glass composition and environment, making the comparison of different 
experiments more subtle, especially if we take into consideration the difference between 
the behaviour in water and in moist atmosphere. On the other hand, several conclusions 
were based on the indirect effect of these phenomena on the measurements of strength 
of glass samples subjected to different kinds of treatments (Han and Tomozawa 1989). 
A condition which has often been shown to lead to several superposed effects on the 
flaw distribution and nature in the glass, but also on the local state of the residual 
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uniform dissolution 




Figure 11. Schematic representation of the blunting mechanisms (from Ito and 
Tomozawa 1982). 

stresses which strongly affects these types of measurements. For example, in the case 
of the aging of indentation flaws, the experiments of Lawn et al (1985) suggest that the 
strength increase during aging is caused by the relaxation of the indentation residual 
strains (also produced to a variable extent when finishing the glass surface) through 
subcritical crack propagation of the flaws. 

Fett et al (2005b) have recently modelled quantitatively the strength increase 
of an uncharged soda-lime glass plate after inducing a surface compression layer by 
ion exchange in water at 90 C, and obtained good agreement with experimental 
measurements. They also made an analogue model for the effects of the formation 
of an ion exchange layer in the internal crack walls and tip neighbourhood of both static 
and propagating cracks in a soda-lime glass in water (Fett et al 2005a). This model 
supports the interpretation of the threshold behaviour by the effect of the compressive 
stresses induced in the ion exchanged layer. They finally compared (by AFM post- 
mortem recombination) the mismatch in the morphology of the crack surfaces obtained 
after inducing the partial corrosion of the crack walls of an arrested crack in soda- 
hme glass (Guin and Wiederhorn 2005) and letting the crack repropagate to failure. 
The comparison reveals a homogeneous corrosion of the two faces and is argued to be 
evidence for the sharp nature of the crack tip at the propagation threshold. 

Although these measurements are quite convincing for the specific case of soda-lime 
glass, the applicability to all conditions requires more systematic work of this kind. It 
is generally accepted that after extreme corrosive attack (such as after HF etching) the 
cracks are very likely to be blunted as proposed by Proctor (1962). However, several 
intermediate conditions may exist where the elevated susceptibility of some kinds of 
glass to dissolution induces an anomalous propagation behaviour as in the high alkali 
binary glasses studied by Gehrke et al (1991) or in some alumino-phosphate glasses 
(Etter et al 2008). 
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4.3. Crack tip plasticity 

Another major objection to the sharp-crack atomic-bonding paradigm comes from the 
hypothesis of a significant contribution of plastic deformation at the crack tip in the 
slow crack propagation mechanisms. Since the early works of Dugdale (1960) on metals, 
and the observations of plastic behaviour of glass under compression in indentation and 
scratch marks (Taylor 1949), several investigations have been carried out to track the 
evidence of plastic behaviour in the strong tensile stress field at crack tips in glasses. In 
particular, Marsh (1964a; 1964b and refs there in) edited a series of papers developing a 
complete theory of elasto-plastic fracture propagation for glass, which was proposed as 
an alternative to the Griffith explanation for brittle glass strength and to the Charles- 
Hilligs interpretation of static fatigue as the effect of stress-corrosion. 

The theory closely follows the formahsm developed by Dugdale for metals, defining 
as a suitable fracture criterion the formation of a plastic zone of critical dimensions ~ 6 
nm (Marsh 1964b) for typical glasses having an estimated yield stress of about 6 GPa. 
The occurrence of static fatigue is explained as a contribution of the stress-enhanced 
water diffusion into the glass to a local reduction of the glass viscosity and thus of the 
resistance to plastic ffow (Marsh 1964b). The time scale for fracture propagation is thus 
attributed to the transport of water in bulk glass and to the consequent viscous fiow, 
although the kinetics is not explicitly developed. 

Wiederhorn (1969) performed accurate measurements of the fracture energy for six 
glasses in inert environment and showed that once the stress-corrosion is suppressed, 
the fracture surface energies 7/ = Gc/2 range between 3.5 and 5.3 J/m^, a figure which 
is ten fold higher than the typical values of the surface tension of glasses 70 ~ 0.5 J/m^ 
estimated by Griffith (1921) and which suggests a major contribution of irreversible 
processes in crack propagation. By analysing the implications of the plastic model 
proposed by Marsh (1964b), Wiederhorn (1969) estimated a plastic work of 4.5 J/m^ 
for a plastic zone of 2.6 nm and a yield stress of 10 GPa, and he observed that although 
it is in rough agreement with the measured fracture energies it is very small compared 
to that obtained in metals of 10^ J/m^. Moreover, since all the estimated space scales 
for the local deformations are close to the dimensions of the molecular scale (the average 
size of glass network rings is of the order of 0.5 nm), he questioned the opportunity of 
applying the concept of 'plasticity' to the nanoscale. 

On the other hand, the sharp-crack atomic-bonding paradigm that was advocated 
by Lawn (1983) as the foundation of the brittle fracture theory, considers that load is 
borne in a completely elastic way at the molecular level. According to this hypothesis 
only the molecules at the direct end of the molecularly sharp crack tip are broken 
sequentially under the elevated load. Lawn et al (1980) conducted an extensive 
evaluation of this possibility by comparing the nature of the atomic bonds and structure 
between glasses and other brittle or ductile materials. They supported their conclusion 
on the absence of plastic behaviour in glasses demonstrating by TEM observations the 
absence of dislocations in crystalline oxides with same composition (as opposed to more 
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ductile crystals). 

However, the hypothesis of Marsh (1964b) went well beyond the discussion of 
fracture in inert atmosphere and vacuum conditions. Marsh was well aware that water 
is the oxide which has the highest effect in lowering the viscosity of a glass melt (cf. 
Del Gaudio et al 2007) and his explanation of the effect of humidity on the slow crack 
propagation (as an effect of reduced glass viscosity by stress-enhanced diffusion of water 
into the glass) was intended as a completely alternative mechanism for stress-corrosion. 
By looking at this interpretation with more modern eyes, the basic mechanism remains 
indeed similar, since the reduction of glass viscosity is probably caused by a stress- 
enhanced hydrolysis of the network siloxane bridges, substituted with more mobile 
SiOH terminations. However, this hydrolysis would not be limited to the crack tip 
rings, but rather involve a larger process zone, of the order of a few tens of nanometres, 
corresponding to a region with an effective yield stress of 1 GPa. The model is also 
quite different in that the time scale is not simply governed by the rate of reaction at 
the tip, but also by the rate of stress-enhanced diffusion of water into the glass and the 
time scale of structural relaxation which is accelerated by both stress and the reduced 
viscosity due to hydration. This would lead to a process zone size that is a decreasing 
function of the crack velocity and thus of K. 

The development and experimental investigation of a similar scenario was 
progressively carried on by Tomozawa and collaborators in the eighties using IR 
and Raman spectroscopies to detect water diffusion profiles in silica glass and the 
consequent changes in the glass structure, identified by changes of the so called 'fictive 
temperature Q- A good review on this work can be found in Tomozawa (1996). However, 
some relevant points will be recalled here. Ito and Tomozawa (1981) initially observed 
an increase in the dissolution rate of silica glass under hydrostatic pressure at 285°C, and 
argued that the opposite should happen in tension, in contrast to the stress-corrosion 
hypothesis of Charles and Hilligs (1962). On the other hand, Nogami and Tomozawa 
(1984) showed that the water diffusion into glass was exponentially enhanced on the 
tensile side of a bent plate of silica glass under saturated moist atmosphere at 192 °C 
(the opposite occurring on the compressive side), the hydrated layer reaching a thickness 
of several microns after 300 h. These two combined observations led Tomozawa to 
suggest that glass structural relaxation caused by stress-induced water penetration at 
the crack tip could be an alternative mechanism for slow crack propagation in glass. 

+ Glass being an out of equilibrium material, its structure and properties are not uniquely determined 
by the thermodynamic variables, but also depend on the past thermal history. According to the 
hypothesis of Tool (1946), the thermal history may be represented by the knowledge of a unique 
additional variable named the fictive temperature Tf. This variable represents the temperature of 
the supercooled equilibrated glass melt which has the same structure as the present glass. For the 
same composition, a higher fictive temperature corresponds in general to a lower density and index of 
refraction. When a glass is aged at a temperature T' < Tf, its structure undergoes a slow evolution 
toward the equilibrated state, the relaxation time being inversely proportional to the glass viscosity. 
During this relaxation the fictive temperature decreases towards the value of T'. However, in practice 
it is quite difficult to equilibrate glasses to a fictive temperature lower than 900 °C. 
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However, most measurements were done at high temperature and the extrapolation to 
ambient condition is generally difficult due to non equilibrated reactive diffusion and 
to difficulties in estimating the solubility, especially due to reported trend inversions at 
different temperatures. A major advance was provided by Tomozawa et al (1991) who 
measured by resonant nuclear reaction analysis the enhanced water entry into silica glass 
fracture surfaces during slow crack growth at ambient temperature. The thickness of the 
hydrated layer at ambient condition was estimated at about 10 nm which is of the order 
of the experimental resolution, but the increase in the amount of water is significant 
when compared to crack propagation in inert environment. The support for the overall 
model has progressively advanced through several complementary observations, which 
are necessary to separate and test the effect of the water content, the fictive temperature 
and the applied stress. Peng et al (1997) have documented the effect of stress and water 
content in accelerating the surface structural relaxation in silica fibres. Li et al (1995) 
have shown that silicate glasses with higher fictive temperature present both greater 
fatigue resistance and greater inert strength. Koike and Tomozawa (2006) showed that 
the fictive temperature also affects differently the subcritical crack growth rate of silicate 
glasses and soda-lime glasses. 

The bibliography on this subject is very rich and can not be reviewed in detail 
here. Several observations and especially some interpretations have been long debated 
and are still disregarded by the 'sharp-tip' community, since once again they are based 
on many indirect observations and extrapolations from different temperatures and other 
conditions. Moreover, the thermal treatments to change the fictive temperature are not 
without effects on the flaw distributions and the effects on strength should be interpreted 
with great care. However, the Tomozawa mechanism remains largely plausible and 
would need a major synthetic effort to prove its real effectiveness for crack propagation 
in silica (and especially in other glasses) at ambient condition. Such modelling should 
end in the writing of kinetic equations that would fit the experimental data as well as 
the Wiederhorn equation (fl]) and that would provide a physical interpretation for the 
model coefficients related to measurable quantities. 

A different series of investigations were made possible by the application AFM 
techniques in glass science in the '90s. In 1996 Guilloteau et al made the first in-situ 
AFM measurements of the external glass surface during crack propagation in borosilicate 
glass. They observed the presence of a surface depression ahead of the crack tip, and 
they estimated the presence of a process zone of 50 nm size (corresponding to a yield 
stress of 1 GPa) by identifying the deviation of the vertical displacement profiles from 
the power law predictions of linear elasticity. Similar results were obtained by Celarie 
et al (2003) on lithium alumino-silicate glasses and by Prades et al (2005) on silica 
glass, by in-situ AFM observation of slow fracture propagation in DCDC samples under 
pure mode I. The process zone size was observed to grow in size from 20 to 100 nm 
when the propagation velocity diminished from 10~^° to 10~^^ m/s due to a decrease 
of K. Moreover, these AFM measurements suggested that the crack propagation in the 
process zone proceeds by the nucleation, growth and coalescence of nanometric cavities. 
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in a similar way to what happens in the ductile fracture of metals at the micrometer 
scale. 

These observations reopened the debate on crack tip plasticity, which is at present 
very active due to several contradictory observations. Lopez-Cepero et al (2007) 
have observed that the post-mortem recombination of the AFM measurements of the 
morphology of opposing crack surfaces in silica and soda-lime showed no evidence for 
the expected traces of the nanocavities in the bulk of the specimen. Moreover, Fett 
et al (2008) showed the inadequacy of the 2D plane-stress solution of Maugis (1992) 
used by Guilloteau et al (1996) and then by Celarie et al (2003) to represent the elastic 
vertical surface displacement near the crack tip with a r~^/^ dependence. When using 
the correct r'^ dependence (with A ^ 1/2), which was derived by Bazant and Essentoro 
(1979) to describe the corner singularity of a 3D surface striking crack (cf. Dimitrov et 
al 2006), the deviation of the displacement profiles measured by AFM from the elastic 
solution is much less evident (Fett et al 2008). The quality of the fit on the vertical 
topography profiles is deeply affected by the presence of the residual roughness even 
after excellent polishing. Moreover, the combination of the roughness and the crack 
tip surface deformation could produce very similar features to the cavities observed by 
Celarie et al 2003 as shown by numerical simulation (Fett et al 2008). 

However, all these measurements are very close to the instrumental limits of the 
AFM topographical images which are made with AFM tips of an average radius of 
10 nm. Better insights could in principle be obtained through the delicate approach 
of separating the vertical surface displacement from the surface roughness by using 
digital image correlation techniques (Hild and Roux 2006). Moreover, a signature of the 
occurrence of inelastic processes during crack propagation can be identified through the 
statistical correlation functions of the crack surface morphology (cf. section 

The recent development of complementary techniques to measure the crack tip 
stress fields with nanometric resolution are also very promising for determining the 
effective scale of a damage process zone. We can cite the SEM cathodo-luminescence 
measurements (Zhu et al 2007; Pezzotti et al 2008) and the promising measurements 
by nanoRaman and electron back scattered diffraction techniques (Vaudin et al 2008). 

4.4- Local crack tip environment 

Since crack propagation in the stress-corrosion regime is strongly affected by the 
environment at the crack tip, and since this is a highly confined region with a typical 
crack opening ranging from 0.5 nm at the tip to a few nm at a distance of 1 /xm, it is 
questionable whether the local crack tip environment can be treated as equivalent to 
the outer environment, be it liquid or gas, and also whether the macroscopic description 
of the thermodynamic state of the environment at a nanometric confined crack tip is 
relevant. 

The region II behaviour is a first sign of an inhomogeneous crack tip environment. 
Since water is being used by the crack tip for propagation, the local concentration (in 



CONTENTS 



25 



gas or solute) at the crack tip is considered to decrease when the crack tip velocity is 
comparable with the transport rate of water in the cavity. The near tip concentration 
in region II was modelled to fall to near zero (Wiederhorn 1967), and the concentration 
gradient established in the crack cavity would be the motor for water diffusion to the 
crack tip, its rate becoming the controlling factor for crack propagation (cf. section 
I3.2p . On the other hand, when the whole sample is plunged into water, a more subtle 
issue is given by the viscous transport of liquid water toward the tip. This transport is 
indeed very effective since the region I then only starts changing its slope when near to 
region III. However, at elevated velocity, a pressure drop at the tip may be caused by the 
increasing viscous drag, leading eventually to a negative pressure in the crack tip region. 
A very interesting investigation by Michalske and Frechette (1980) has suggested that 
a strongly negative pressure can be sustained at the crack tip without cavitation due to 
the extreme confinement, leading to the observation of a supercritical crack propagation 
at reduced velocity, which can be interpreted according to the shielding effect (reduction 
of K*) caused by the internal forces in the crack cavity (cf. Fig. [1]). A further increase 
of K eventually leads to cavitation in water when the negative pressure extends to a 
less confined region of the crack cavity, followed by an instantaneous order of magnitude 
rise in the crack velocity due to the sudden depletion of the crack tip environment and 
to the release of the internal forces. A similar cavitation-induced dynamical instability 
was observed by Maugis (1985) in brittle polymers. 

In alkali containing glasses in water, the enhanced leaching at the crack tip may 
cause a local change of the hquid composition at the crack tip. This was suggested to 
be the cause of an observed transition in the slope of the region I in soda-lime glass 
(Wiederhorn and Johnson 1973). At low crack velocity the crack tip environment would 
be equilibrated to the external liquid composition, while becoming more related to the 
crack tip reactions at higher velocity. The same kind of effect can be supposed for other 
type of corrosion products, such as silicic acid, but still has to be quantified. 

Recent AFM in-situ observations have shown that nanometre scale capillary 
condensation occurs at crack tips in silica glass in moist air for crack velocities between 
10"^2 and 10"^ m/s (Wondraczek et al 2006; Ciccotti et al 2008). This hquid phase 
is made stable by the short range interaction with the glass surfaces under elevated 
confinement and constitutes a major alteration of the local environmental condition 
at the crack tip with respect to the gaseous hypothesis (Wiederhorn 1967). In this 
condition, not only is the water readily available at the crack tip, but the local pressure 
becomes strongly negative due to the Laplace pressure, thus exerting a strong attractive 
force between the crack lips and altering the equilibrium of the stress-corrosion reactions. 
Since the crack-tip reaction rate is expected to depend on the chemical activity of 
water, and the condensate was shown to be in equilibrium with the humid atmosphere 
(Grimaldi et al 2008), we can expect that the humidity dependence of the crack 
velocity remains similar to what expected in gaseous environment, yet the situation 
should change when transport phenomena become relevant or when the equilibrium is 
questioned. For example, the liquid condensate can be an explanation of the reduction 
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of the transport limited region II region observed at elevated humidity (Wan, Aimard et 
al 1990). The limiting factor for the diffusion of water is shifted from the very confined 
crack tip region to the end of the liquid condensate where the larger crack opening allows 
faster diffusion. 

The liquid condensation was also observed on alkali containing glasses (Celarie et 
al 2007) and should in this case have a much larger effect since the local compositional 
changes, due to the crack tip enhanced leaching, can not be diluted like in a bulk 
liquid. Preliminary results show that the extent of the condensation in soda-lime glass 
is significantly enhanced by the alterations due to the local leaching (cf. section H75|) . 
The relevance of capillary condensation in the propagation kinetics of phosphate glasses, 
which are less resistant against dissolution, can be seen in the anomalous temperature 
behaviour observed in the v{K) relations measured by Crichton et al (1999) and similarly 
by Etter et al (2008) at ambient temperature. The anomaly is suppressed at higher 
temperature where the condensation can not form. 

4-5. Alkali ion migration under stress gradient 

In section 14.11 we have seen that chemically driven sodium ion migration can happen 
spontaneously at soda-lime glass surfaces thus inducing a local state of compression in a 
surface hydrated layer. Gorsky (1935) first showed that the opposite can also be true, i.e. 
that a strong stress gradient can induce sodium migration towards the tensile direction of 
the gradient. While the spontaneous slow flow is generally balanced by interchange with 
other positive ions (such as the hydronium ion) in order to preserve charge neutrality 
(Doremus 1975; Lanford et al 1979), the behaviour under a strong stress gradient is less 
evident. Weber and Goldstein (1964) measured a transient electric current between the 
tension and compression sides of a bent soda-lime slide, indicating that the non-balanced 
flow can occur before being counteracted by the build up of an opposite electric potential 
difference. Moreover, Langford et al (1991) have revealed intense sodium emission after 
rapidly fracturing soda-lime glass in vacuum. The delay in the emission was of the 
order of a few tens of milliseconds, indicating that sodium migration can be very fast 
in the presence of strong stress gradients and can thus also affect region III water-free 
fast crack propagation. Several observations were made concerning an excess sodium 
concentration in the outermost surface layer of fracture surfaces in so da- lime, overlaying 
a sodium depleted region of larger thickness of hundreds of nanometres (Pantano, 1985). 
This phenomenon looks similar to the chemically driven migration, however its space 
and time scales are quite different and more investigation is needed concerning the 
interaction of this phenomenon with the external environment during and after the 
fracture process. 

Watanabe et al (1994) made an interesting AFM investigation of the time evolution 
of the fracture surfaces of a soda-lime glass broken in low vacuum at a velocity of 10~^ 
m/s (region III) and then aged in ambient moist atmosphere. The excess concentration 
of sodium on the fracture surfaces is the cause of an unusually rapid corrosive action on 
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these surfaces, observable by the appearance of protuberances and swelhngs interpreted 
as the recondensation of corrosion products in the form of a weak sodo-sihcate gel 
material. The evolution of such swellings under a common atmosphere finally evolves 
to the formation of carbonatic crystallites due to the interactions with CO2 dissolved 
in the sodium rich wet layers. 

A similar kind of alteration has been observed by imaging external surfaces in the 
neighbourhood of a crack tip produced by indentation on a soda-lime slide (Nghiem 
1998). More recently, Celarie et al (2007) have conducted a systematic in situ AFM 
study of the kinetics of the growth and evolution of the protuberances during crack 
propagation in soda-lime glass under controlled atmosphere of nitrogen and variable 
humidity levels. The absence of CO2 allowed focusing on the space and time scales 
of the sodium exchange mechanisms and of the surface corrosion process. The region 
affected by the surface swellings is shown to have a parabolic shape which is modelled 
according to a competition between the spreading of the surface reaction and the crack 
propagation velocity (Fig. [T2l) . The parabolic shapes for different crack velocities at 
45% RH are consistent with a diffusive process with an effective diffusion coefficient 
of the order of 1 nm^/s, which is an increasing function of the relative humidity. The 
application of phase imaging techniques has also allowed to identify the thickening 
of the sodium enriched liquid layer (predicted by Watanabe et al 1994) prior to the 
manifestation of the swellings. 
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Figure 12. AFM topographical images of a propagating crack in soda-lime glass at 
different crack velocities for RH=45%. 

When fracture velocity slows down below 10^^ m/s the ion exchange front can 
spread ahead of the crack tip and thus start relaxing the tensional stresses in the crack 
tip region, which in turn will cause the progressive arrest of the crack propagation. This 
suggests the extension to the propagation in moist atmosphere of the interpretation of 
the threshold behaviour proposed by Michalske and Bunker (1985). The enhanced liquid 
condensation stimulated by the sodium enrichment at the fracture surfaces and in the 
parabolic front plays the role of the bulk liquid reservoir, but the strong confinement is 
likely to cause different behaviours, like in the observations of Watanabe et al (1994). 
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The space and time scales for the spreading of the ion exchange can be quite different. 

Fett et al (2005a) modelled the mechanical effect caused by the chemical induced 
ion exchange layer forming on the crack surfaces and at the crack tip of arrested and 
propagating cracks in soda-lime glass. This modelling is based on the ion exchange 
kinetics proposed by Lanford et al 1979 for an etched soda-lime surface in water at 
pH7. The lateral extension of the parabolic front predicted by this model would be 
about 1000 times smaller than the one observed in Fig. [121 The extremely different 
scales can be attributed on one side to the discussed difference of the propagation in 
moist air (and the consequent change in the composition of the condensed liquid layers), 
and on the other side to the influence of the strong stress gradient at the crack tip, which 
is not accounted in the model of Fett et al (2005a). 

4-6. Fracture surfaces 

As we have seen throughout this review, post-mortem study of the fracture surfaces 
is an excellent complementary tool to understand the physics and chemistry of the 
mechanisms of crack propagation. On one side, the morphology of the fracture surface 
can contain relevant information on the dynamics of crack propagation and on the traces 
of inelastic processes or of the corrosion at the crack surfaces. On the other side, modern 
high resolution chemical probes can provide invaluable insights into the nature of these 
high energy surfaces and of the mechanisms that originated them (Pantano 1985). 

The fracture surface morphology has long been used to draw information about 
the origin and propagation of the cracks that led a glass specimen to failure. From 
an optical inspection, one can distinguish three typical regions of increasing roughness 
called 'mirror', 'mist' and 'hackle' which correspond to progressive crack acceleration 
from subcritical to dynamic propagation (Johnson and Holloway 1966). The spread of 
AFM investigations in the '90s has shown that surface roughness in the 'mirror' region 
is of the order of a fraction of nm (RMS) when measured on micrometer size images, 
and that it is almost comparable with the roughness of melt surfaces and with the best 
levels of surface finishing by polish treatments (cf. Radlein and Frischat 1997; Arribart 
and Abriou 2000). 

The dependence of the RMS roughness of fracture surfaces in the mirror region 
on the nature of glass was extensively studied to understand the relation between the 
heterogeneities in the glass composition and the crack path through the glass (Gupta 
et al 2000; Wiederhorn et al 2007). Other investigators have explored the scaling 
properties of the fracture surfaces (Bouchaud 1997) and the dependence of the self- 
affine exponents and the characteristic cutoff length on the velocity of crack propagation 
(Bonamy et al 2006). 

The extensive development of dynamical models for the crack front line in a 
heterogeneous material (cf. Ramanathan et al 1997), have made the investigation of 
the scaling behaviour of glass fracture surfaces a very interesting tool for probing the 
physics of crack propagation. A particularly interesting aspect is the correlation between 
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the position of some cutoff lengths in the scahng properties with the size of the process 
zone as observed in several different materials (Bonamy et al 2006). These observations 
support the idea that a generalised process zone behaviour can affect scales of the order 
of a few tens of nanometres during slow crack propagation in silica glass, the size of the 
process zone being an increasing function of the decreasing crack propagation velocity. 

A useful complementary tool to study the damage mechanisms induced by crack 
propagation is the analysis of the mismatches in the recombination of the two opposite 
fracture surfaces. Guin et al 2004 showed that it is possible to match down to nanometre 
scale opposite fracture faces of a soda-lime glass broken in water. The fine comparison 
on a few sections normal to the crack surface showed no evidence of damage cavities in 
the bulk of the glass sample (cf. also Fett et al 2008). However, the determination of 
the remnants of damage processes on post-mortem (unstressed) fracture surfaces in glass 
requires extreme care since the relevant information is at the nanometre scale where the 
metrologic capacities of AFM imaging are affected by several delicate problems such 
as drifts, feedback effects, noise and tip shape artifacts (Flemming et al 2005). The 
accurate treatment of these kinds of problems will benefit from the present development 
of digital image correlation techniques (Hild and Roux 2006) and requires the application 
of advanced statistical estimation tools. 

As discussed in section 13.41 post-mortem AFM investigations and recombination 
techniques have also led to important insights into the nature of the propagation 
threshold and the mechanisms of crack arrest and aging. For example, Guin et al 
2005 could measure the extent of the dissolution of the crack surfaces when a crack in 
soda-lime glass is arrested and aged near the propagation threshold in a basic solution. 
The work of Watanabe et al 1994 on the evolution of the fracture surfaces of a soda-lime 
glass after fracture in vacuum has provided important insight into the time scales of the 
corrosion mechanisms caused by the interaction with the moist atmosphere. 

The UHV AFM investigations of the fracture surfaces of different kinds of glasses 
with molecular resolution (Frischat et al 2004) have provided an intimate view of the 
complex nanostructures in multi-component glasses and their role in the slow crack 
propagation processes. These kinds of measurements are very promising for relating the 
structural measurements in bulk glasses (cf. Wright 2008) with the molecular simulations 
of the peculiar structure of the glass fracture surfaces (Du and Cormack 2005; Roder et 
al 2001). 

Many techniques for the high resolution characterisation of the structural and 
chemical properties of glass surfaces are being developed, which are out of the scope of 
the present review. I will just cite here some relevant investigation of the hydroxylation 
behaviour of fracture surfaces (Souza and Pantano 2002) since the study of the water 
adsorption sites and of the time scales of the hydroxylation mechanisms are of particular 
interest for understanding the high energy of the glass fracture surfaces (Leed and 
Pantano 2003) and modelling their wetting properties, which were shown to have a 
paramount role in the corrosion mechanisms in moist environment (cf. section 22]) • 
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5. Perspectives 

After almost one hundred years from the initial insights into the nature of brittle crack 
propagation in glasses we finally feel the flavour of a turning point in the comprehension 
of the basic mechanisms of slow crack propagation in glasses. The kernel of the difficulty 
turned out to be that there are many actors playing in the opera, most of them acting 
at the nanoscale in a deeply interrelated manner. The amorphous nature of glass has 
deprived the scientists of most of their typical modelling and characterisation tools, 
and the exceptional homogeneity of glass relegated all the relevant phenomena to the 
nanoscale, where the chemical nature of surfaces becomes most influent and the strong 
conflnement makes the macroscopic thermodynamics fail. On the other hand, the 
disordered structure of both the glass network and the reactive molecules (such as water) 
tend to appear ordered at the nanoscale, a fascinating order that is often speciflc to each 
of the many glass compositions. 

The turning point stems from the progressive adaptation to the nanoscale of 
both the investigation tools and the simulation techniques, along with the increasing 
comprehension of the physics of the nanoworld. Different kinds of local scanning probes 
have allowed the visualisation of the space and time scales of the physico-chemical 
mechanisms during or after crack propagation. Several structural and compositional 
analyses arc also increasing their resolution and reducing the severity of their condition 
of use. The most promising developments concern the combination of different 
complementary probes on a single setup (such as AFM+micro Raman techniques). At 
the same time, the gap between simulation and experience is also being drastically 
reduced. On one hand, both molecular dynamics and flrst principle simulations permit 
the modelling of increasing volumes of matter and increasing time durations (which 
remain however the limiting factor). On the other hand, several experiments under 
UHV conditions are exploring physical conditions that approach those being simulated 
(by exploring for example the interaction of few water molecules with the glass fracture 
surface) . 

However, when trying to model these new exciting investigations, we should not 
forget all the ingredients of a sound mechanical modeUing and the subtle details of glass 
surface chemistry. The crack tip and process zone phenomena should be adequately 
separated from the structural behaviour of the glass sample in order to tackle correctly 
the issues related to the management of residual stresses, to surface flnishing, and 
to the effects of the conditioning treatments on the population of surface flaws and 
nucleation sites. A major difficulty highlighted in this review is the subtle difference 
between the properties of crack propagation or surface corrosion processes in liquid 
solutions and under moist atmosphere. While the first case corresponds to an overall 
controlled environment (at least for slow crack propagation) , in the second case the local 
environment at the crack tip or at the glass fracture and external surfaces is determined 
by the evolution in time of a very confined liquid condensation. Carbon dioxide in the 
atmosphere was also shown to be a major actor in the evolution of the different chemical 
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reactions. It is thus of fundamental importance to dispose of a high degree of control 
of the atmosphere during the investigations, and to repeat the experiences under pure 
mixtures of nitrogen and humidity. 

This review is not exhaustive, but hopefully sufficient analysis has been provided 
in order to promote a deep reflection on the relevance of the different actors. This 
reflection should guide to the identiflcation of the space and time scales relevant for 
the action of the different physical and chemical processes in the close neighbourhood 
of the crack tip and crack surfaces as a function of the different glass compositions 
and controlled environments. This analysis should always be complemented by the 
estimation of the energetic contribution of each individual mechanism. Special attention 
should be paid to the systematic complete characterisation of the glass samples used, 
including their detailed composition, structural and mechanical properties and to the 
definition of the thermal history as well as the details of the surface finishing techniques. 
The reproducibility of the results in different laboratories will be key point for a real 
understanding. 

We should not forget that understanding the nature of the crack propagation 
mechanisms will help in solving just part of the problem of glass strength. This line of 
research should be combined with a comprehension of the flaw nucleation mechanisms 
starting from different kinds of defects within the glass and with the development of 
techniques making glass objects more tolerant to the presence of such flaws. Strong 
efforts in these directions are being invested by both the academic and industrial 
communities in order to make glass a promising material for advanced applications 
in the XXI century. 
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